Skip to content

Conversation

@sbrus89
Copy link

@sbrus89 sbrus89 commented Apr 22, 2025

This PR adds the design document for the pressure gradient module for Omega-1. Compiled version here.

Checklist

  • Documentation:
    • Design document has been generated and added to the docs
    • Documentation has been built locally and changes look as expected

Copy link
Collaborator

@vanroekel vanroekel left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good, just a couple tiny questions from me.

Copy link

@philipwjones philipwjones left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Just a couple of questions...

@mark-petersen
Copy link

Test compiled and posted here.

Copy link

@xylar xylar left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@sbrus89, a few things to chat about at our Wednesday meeting (or sooner if you prefer). Nothing that should hold up moving forward on this.

## Verification and Testing

### Test: Spatial convergence to exact solution
For a given analytical $v$, $h$, and $b$, the spatial convergence of the pressure gradient can be assessed by computing errors on progressively finer meshes.
Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What is $b$ in this case?

Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think maybe this should be fleshed out a bit more. What analytic test would we use?

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I revised this part slightly. I think we just need to specify arbitrary functions for the $\alpha$, $h$, and $z$, then we can test the convergence of the discrete tendency compared to the analytical solution.

@sbrus89 sbrus89 force-pushed the omega/pgrad-design branch 4 times, most recently from 674c345 to 16c3c49 Compare September 17, 2025 14:05
@mark-petersen
Copy link

Latest version compiled here.

Copy link

@mark-petersen mark-petersen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for your work on this. It looks great!


### 2.2 Requirement: Initial support for a simple centered pressure gradient

For initial global cases without ice shelf cavities, the pressure and geopotential gradients will be computed with a simple centered difference approximation. In later versions of Omega, one or more high-order pressure gradients will be implemented and will replace the centered approach in production runs.
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

one or more seems different from the next requirement where Adcroft 2008 is mentioned. Is the plan to leave open for multiple options?

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

the hedge Xylar proposed below seems like a good choice will->may and resolves the inconsistency

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Projects

None yet

Development

Successfully merging this pull request may close these issues.

7 participants